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r^ • We discuss the low-temperature dynamics of magnetic flux lines in high-temperature 

fj I superconductors in the presence of a family of parallel twin planes that contain the c axis. 

A current applied along the twin planes drives flux motion in the direction transverse to 

^ ■ the planes and acts like an electric fleld applied to one- dimensional carriers in disordered 

d • semiconductors. As in flux arrays with columnar pins, there is a regime where the dynamics 

is dominated by superkink excitations that correspond to Mott variable range hopping 

(VRH) of carriers. In one dimension, however, rare events, such as large regions void of 

twin planes, can impede VRH and dominate transport in samples that are sufflciently 

long in the direction of flux motion. In short samples rare regions can be responsible 

for mesoscospic effects. The phase boundaries separating various transport regimes are 

discussed. The effect of tilting the applied field out of the twin planes is also considered. 

In this case the resistivity from flux motion is found to depend strongly on the tilt angle. 
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1. Introduction 

The static and dynamical properties of magnetic flux lines in copper-oxide supercon- 
ductors are strongly affected by pinning by point, linear and planar disorder [|l|. Twin 
boundaries are an example of planar disorder that is generally present in superconducting 
YBa2CusOr-x and La2Cu04, where they are needed to accomodate the strains produced 
by a crystallograflc tetragonal-to-orthorhombic transition. Twins most often occurs in two 
orthogonal families of lamellae forming a mosaic . It is also possible to prepare samples 
that contain a single family of parallel twin planes |]^,Q . Columnar or linear defects can be 
produced by bombardment of the crystal with energetic heavy ions or are embodied in 
forests of screw dislocations parallel to the crystal growth direction . Both twins 0,^§] 
and columnar defects Q constitute examples of macroscopic correlated disorder that can 
be responsible for a sharp decrease in the resistivity for specific field orientations. 

Extensive investigations of twin-boundary pinning have been carried out by Kwok and 
coworkers [^. These authors studied a variety of YBCO single crystal samples containing 
single families of parallel twins lying in planes spanned by the c axis, with spacings ranging 
from microns down to several hundred Angstroms. As indicated by earlier decoration 
experiments @ , the order parameter is suppressed at a twin boundary at low temperatures 
and the twin attracts or pins the vortices. Transport experiments show clear evidence of 
strong twin-boundary pinning even in the flux liquid phase 0. The pinning is most 
effective when the external fleld is applied along the c axis and the driving current lies in 
the ab plane and is parallel to the plane of the twins, resulting in a Lorentz force normal 
to the twin planes (Fig. la). In this case twin planes act like strong pinning centers 
and the linear resistivity as a function of temperature drops sharply at a characteristic 
temperature that marks the onset of twin-boundary pinning . In samples with only a few 
widely spaced twin planes this drop in the resistivity is followed by an abrupt shoulder at 
a lower temperature, corresponding to the flrst order freezing transition into an Abrikosov 
lattice. The abrupt shoulder is not observed in heavily twinned samples where the freezing 
transition is apparently suppressed by disorder. Additional experimental evidence for twin- 
boundary pinning comes from the observation of a sharp downward dip in the resistivity 
as a function of angle as the external fleld is rotated out of the twin plane through the 
c direction Q (see Fig. lb). This strong angular dependence is a clear signature of 
anisotropic pinning by twin planes since point disorder in the form of oxygen vacancies, 
while certainly present in all these samples, can only yield a weak dependence of the 
resistivity on the tilting angle. 
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Another experimental probe of flux-line dynamics that has been recently used to study 
the effect of correlated disorder is real-time imaging of flux proflles |jiy,|ll|. Imaging ex- 
periments of the fleld penetrating into single crystals with families of twins lying in planes 
spanned by the c axis (for flelds directed along the c axis) show strong pinning by twin 
planes for flux motion in the direction normal to the planes, conflrming the strong twin- 
boundary pinning seen by transport measurements for Lorentz forces applied normal to 
the twin planes. Imaging of flux motion along the twin planes have, however, given con- 



tradictory results. Duran and coworkers [|T^ argued that flux penetrates more easily in 



the twin regions than in the channels between the twin planes, in apparent contrast with 
the observation of a lowered flux-flow resistivity by twin-boundary pinning by the Argonne 
group for the parallel geometry, where the Lorentz force is applied along the plane of the 
twins [0. A more recent imaging experiment by Vlasko-Vlasov and collaborators [Tl\ 



shows, however, that the twin planes act as planar pinning barriers for all directions of 
flux flow, giving rise to guided vortex motion. In this paper we are only concerned with 
transport with flux motion in the direction transverse to the twin planes. For this geom- 
etry all experimental probes of flux dynamics conflrm that twin planes provide attractive 
potential wells for the vortices. A detailed theoretical understanding of the imaging ex- 
periments and their relationship to transport transport in the parallel geometry remains, 
however, an open question. 

The static and dynamical properties of flux line assemblies in the presence of a random 

array of columnar pins have been studied in detail by mapping the physics of magnetic flux 

lines onto the problem of localization of quantum mechanical bosons in two dimensions 

T^ . This mapping exploits methods developed to understand the behavior of ife^ fllms 



on disordered substrates [|IJ] and to decribe electronic transport in disordered supercon- 
ductors Ijl^. At low temperatures there is a "Bose glass" phase, with flux lines localized 
on columnar pins, separated by a phase transition from an entangled flux liquid of delocal- 
ized lines. In the Bose glass phase the linear resistivity vanishes and the current-voltage 
characteristics are nonlinear. Transport in this regime closely resembles the variable-range 
hopping (VRH) of electrons in disordered semiconductors in two dimensions [IB . 



In this paper we employ similar methods to study flux-line dynamics in the presence of 
a single family of parallel twin boundaries lying in planes containing the c axis. For H || c, 
the flux lines are localized by the pinning potential in the direction normal to the planes. 
At low temperatures, when the average vortex spacing ao ~ {(po/B) exceeds the average 
distance d between twin planes, all flux lines are localized on the twins, progressively 



"filling" the planar pins as the field is increased. We only consider vortex arrays that 
are sufficiently dilute (S << S/ ~ (po/cP) that transport is controlled by single- vortex 
creep. This regime is experimentally relevant and has in fact been probed in ffux arrays 



with columnar defects [^ . Flux motion in this regime is dominated by thermally activated 
jumps of the vortices over the relevant pinning energy barriers U{L, J), yielding a resistivity 
p = S/J, given by ^, 

p{T) ^ poe-"^^'-^)/^, (1.1) 

where po is a characteristic flux-flow resistivity. Here we focus on flux motion transverse to 
the twin planes at low fields and temperatures, which resembles the hopping of electrons 
in one- dimensional disordered superconductors. The energy barriers lA corresponding to 
the various low-lying excitations that can contribute to transport are evaluated and are 
summarized in Table 1. Typical phase diagrams in the {J,L) plane displaying the regions 
where the different transport mechanisms dominate are shown in Fig. 2. There is a 
characteristic current scale Jl ~ 1/-^, where L is the sample thickness in the field direction, 
that separates the regions of linear and nonlinear response. As L ^ oo the response is 
always nonlinear. For large enough current the dominant excitations are the half-loop 
configurations shown in Fig. 3a, of transverse width smaller than the average separation 
d between twin planes. For currents below Ji the width of a typical half- loop excitation 
exceeds the mean distance between pins. In this case standard VRH arguments |]TB,ff3 



where an electron hops larger and larger spatial distances to find "good" traps of low energy 
(Fig. 3c), suggest a nonlinear current-voltage characteristic V ~ exp[— (£'fc/T)(Jo/J)^'^] 
in the localized phase, where E^ and Jq are characteristic energy and current scales given 
below. The phase diagrams shown in Fig. 2 are qualitatively similar to the phase diagrams 
one would obtain for vortex arrays in the presence of columnar pins. In that case one finds a 
nonlinear current- voltage characteristic at low currents typical of VRH in two dimensions, 
with V ~ exp[-(Efc/T)(Jo/J)V3]. 

The main difference between ffux motion in the presence of columnar defects and ffux 
motion transverse to an array of parallel twin planes is that in the latter case the low 
temperature transport is one- dimensional. In one dimension VRH can be impeded by the 
presence of large rare regions void of energetically favorable pins [|T^,|T^ . These rare regions 



free of localized states are exponentially rare, but have a very large resistance and can sup- 
press VRH or even dominate transport (Fig. 4) in one dimension since the vortices cannot 
get around them. For samples that are sufficiently long in the direction of fiux-line motion 



so that they contain a large number of such rare regions, this mechanism yields a nonlin- 
ear current-voltage characteristic of the form V ~ (T/E'fc)(J/Jo)^/^ exp[— (£'a;/T)^( Jo/J)]. 
Shorter samples will typically contain only a few rare regions and these will determine the 
sample's resistance. In sufficiently short samples there will be a spread of values of the 
resistance between different samples resulting in reproducible sample-to-sample resistance 
fluctuations. 

Finally, the model presented here is also relevant to the dynamics of Josephson vor- 
tices in artificially structured "giant" Josephson junctions. Consider for instance a planar 
junction of in- plane dimensions large compared to the Josephson penetration length at the 
contact of two superconductors (the plane of the junction is the xz plane) . A magnetic field 
applied in the plane of the junction (say, in the z direction) penetrates into the junction as 
a chain of Josephson vortices which lie in the contact plane [|T9| , |20(| . The intervortex spac- 



ing along the x direction is determined by the strength of the applied field. The vortices 
are localized in the plane of the junction and form therefore a (1 -|- 1) -dimensional vortex 
array in this plane. If the junction is not uniform, the vortices are pinned independently at 
low fields, as indicated by the fact that the critical current does not depend on magnetic 
field. Additional defects can be artificially introduced in the junction. The junction can 
be artificially structured by the introduction of an array of defects along the x direction 
spanning the field axis (z) and the junction thickness (y). At low temperature vortex 
motion along the junction plane will then occur via thermally activated jumps between 
these defects and will be described by the one-dimensional tight-binding model introduced 
below. 

In Section 2 we first review the simple model of interacting fiux lines in the presence 
of correlated disorder introduced by Nelson |2^] and by Nelson and Vinokur [|13|. The 



analogy with quantum mechanics of two-dimensional bosons and the reduction of the low 
temperature dynamics transverse to an array of parallel twin planes to a tight binding 
model in one dimension are then discussed. In Section 3 we estimate the pinning energy 
barriers associated with the low-lying excitations from the ground state and the corre- 
sponding contributions to the resistivity. The phase boundaries separating the regions 
of the (L, J) plane where the various contributions dominate are also discussed. A brief 



summary of these results has been presented elsewhere |]2^ . In section 4 we consider the 
case where the external field is tilted at an angle 9 away from the c axis and out of the 
twin planes. The resistivity displays a strong angular dependence with a sharp downward 
dip at ^ = 0. Finally, in Section 5 we discuss the role of rare fiuctuations in this one 
dimensional geometry. 



2. Vortex Free Energy and Tight-binding Model 

We are interested in transport at low fields and temperatures where each fiux line is 
localized on one or more twin planes. In this regime dominated by single-vortex dynamics 
a detailed description of transport can be developed. Our starting point is a model-free 
energy for fiux lines in a sample of thickness L in the presence of a family of parallel 
twin planes fl^. The field is along the c axis, chosen as the z direction, and the fiux 



lines are parametrized by their trajectories {vi{z)} as they traverse the sample. The twin 
boundaries are parallel to the zx plane. The model free-energy for a single fiux line at 

{vi[z)^ z) is given by 
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with 

M 



Vniy) = J2Vii\y-yk\). (2.2) 

fc=i 

Here Vd is the random potential arising from a set of M x- and ^-independent pinning 

potentials Vi{\y — Yk\) centered at the locations {Yk} of the twin planes. The first term on 

the right hand side of Eq. ( p.l|) is the first term in a small angle expansion of the elastic 

energy of a nearly straight vortex line, with ei ^ {M±/Mz)eo ln{Xab/^ab) the tilt modulus 

and Xab and ^ab the penetration and the coherence lengths in the ab plane, respectively. 

The effective mass ratio M±/Mz << 1 incorporates the material anisotropy and eo ~ 

(^o/47rAa&)^ is a characteristic energy scale. For simplicity we model Vniv) as an array of 

identical one-dimensional square potential wells of depth Uq, width 2bo and average spacing 

d, passing completely through the sample in the x and z directions [^^. Assuming the 

potential wells are centered at uniformly distributed random positions {Yk} and bo « d, 

we find Vd ~ ^0(^)5 while the random potential fiuctuations 5Vd(2/) = Vd(j/) — Vd 

satisfy 

dVD{y)dVD{y')=Ad{y-y'), (2.3) 

with A !^ Uq ^ ^ [l+0{2bo/d)]. The interaction between vortex lines and a twin boundary 
has been studied by Geshkenbein in the context of the Ginzburg-Landau theory p^ . 

In this paper we assume that the twin planes are randomly located along the y direc- 
tion and that the separations between neighboring twins are Poisson-distributed (see also 
Section 5 below) . This appears to be the case in some of the samples employed by the Ar- 
gonne group [|25|. On the other hand, the twin structures that form naturally in YBCO to 



accomodate the strains arising from a tetragonal-to-orthorhombic transformation, which 
take place around lOOOi^ as a result of oxygen vacancy ordering, are often quite different 
0. They consists of lamellae or colonies of parallel twins oriented in either the (110) 
or (110) directions. Orthogonal twin colonies form a mosaic-type structure, containing 
colonies of various size. The colony size scales with the square of the average twin spacing 
d, while the latter remains rather uniform within a given colony 0. There is a repulsive 
interaction between the twin planes of a given colony arising from the stress produced 
by one twin plane in the region of another. This interaction leads to the regular spacing 
of the twins within a colony which resembles an approximately regular one-dimensional 
lattice of twin planes. Vortex dynamics in the presence of such a twin structure will not 
be described by the model presented in this paper. On scales shorter than the typical twin 
colony size vortices are pinned by a regular array of planar defects, while on scales larger 
than the colony size the theory developed to describe vortex dynamics in the presence of 
columnar defects should apply. In contrast, the twin structures observed in the samples 
used by the Argonne group consist of a single colony of parallel twins with large variations 
in the twin spacing. A twin structure of this type may arise if the sample is annealed 
and the twins "fall out of equilibrium" arranging themselves in one-dimensional liquid-like 
structure within a given colony. 

The free energy for an assembly of A^ flux lines is given by 



N „L 

V)=E^^ + 2^/ ^(Ir.(^)-r,(^)|)rf^, (2.4) 



where V{r) is the pair interaction potential, assumed local in z. It can be shown [jT^ that 
both higher order terms in the small angle expansion of the elastic energy and nonlocality 
in z in the pair interaction are negligible provided \dvi/dz\'^ « M^/M^ for the most 
important vortex configurations. In the following we will simply use the form of the pair 
interaction for nearly straight flux lines, 

V[r) ^ 2eo [Ko{r/Kb) - i^o(r/U)] , (2-5) 

where Kq{x) is a modified Bessel function. 

At low temperatures, when the average vortex spacing a^ ~ {(po/B) exceeds the 
average distance between twin planes, all flux lines are localized on the twins, progressively 
"filling" the planar pins as the field is increased. Any real sample will, however, also 
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contain point defects, which are known to promote flux-hne wandering. Recent analytical 
and numerical work has shown that in the case of planar pins a single flux line remains 
localized on the pinning plane even in the presence of additional weak point disorder in 
the bulk, which is always present in real samples [^. The stability of the localized 
Bose glass phase in 1 + 1 dimensions - the case relevant to the model considered here - 
in the presence of point disorder has been studied recently by Hwa et al. ^^■ These 
authors considered a model where flux lines directed along the z direction and confined 
to the zy plane are pinned by the competing action of randomly distributed linear defects 
spanning the plane in the z direction and point defects described by a random potential 
with variance Aq. They showed that in 1 + 1 dimensions the low temperature phase is 
of the Bose glass type with flux lines localized on the linear pins when point disorder 
is weak. The localized phase is marginally unstable to point disorder, but only beyond 
an astronomically large crossover length scale. Point disorder will in general lower the 
energy barriers associated with the various low-lying excitations discussed here. It does 
not, however, have a significant effect on the energy barriers in the rigid flow (or half- loop) 
regime where Urf ~ L (or Uhi ~ ^/J)- This is because the energy gain SF^ associated 
with the pinning of a fluctuation of length L by point defects only grows as L^'^ (or as 
J- 1/2 ij2 the nonlinear regime), with 5Fa ~ (AqL)^'^. Sufficiently strong point disorder 
can, however, lower considerably the barriers for variable range hopping a la Mott since 
both the barrier IAmou and 5F/^ grow like L^'^. In this case to assess whether point or 
correlated disorder dominates one needs to compare quantitatively the relative strengths 
of these two types of crystal defects. This comparison involves unknown parameters and 
is beyond the scope of the present paper. 

The problem of one flux line localized near a single twin plane has been studied by 
Nelson by exploiting the mapping of the statistical mechanics of magnetic flux lines onto 



the quantum mechanics of two-dimensional bosons ||21|]. At zero temperature the twin 
provides a binding energy Uq per unit length for trapping the flux line. In the presence of 
thermal fluctuations Uq is replaced by a smaller binding free energy per unit length U{T), 
to account for the entropy lost by confining the flux line near the twin plane. In infinitely 
thick samples (L ^ oo) this binding free energy is determined by the zero-point energy 
of a fictitious two-dimensional quantum mechanical particle confined to a one-dimensional 
potential well, or U(T) = —Eq{T), where Eq(T) is the ground state eigenvalue of a two- 
dimensional "Schrodinger" equation, 

- i^^l + Viiy) ^o(x, y) = EoM^. y). (2.6) 



and the twin plane is centered at y = 0, i.e., Vi{y) = —Uq for |y| < bo, Vi{y) = for \y\ > bo- 
In the quantum mechanical analogy T plays the role of Planck's constant h, ei that of the 
mass m of the fictitious particle and L~^ that of the particle's temperature. The x and 
y degrees of freedom are decoupled in Eq. (|2.6[ ) and the ground state wavefunction is the 
product of a free-particle wavefunction in the x direction and the ground state wavefunction 
(f)o{y) of a one-dimensional particle in a well [^, ^/jo{x,y) = ^=e*'^^^(^o(2/) [^ , |T3[| . Here 
D is the system size in the x direction and Qx the wavevector. The corresponding ground 
state energy is Eq{T) = -^^ + Eqw The first term is the free particle contribution 
describing the energy cost associated with localizing a fiux line within a distance ~ ^Ti/qx 
and E^yj < is the ground state energy of a one-dimensional particle in a well []2^ . If no 
other fiux lines are present and the sample is infinite in the x direction, we can take 5^=0 
and U{T) = —Eq^ = Uof{T/T*) [^. Here T* = y/2Uoeibo is a crossover temperatiire 
above which thermal fluctuations delocalize the flux line, and f{x) ^ 1 — ^x^ for x « 1 
and /(x) ~ 1/x^ for x » 1. The probability of flnding a point on the vortex at a 
transverse displacement r^ is proportional to |'i/'o(r±)P and depends only on the transverse 
displacement y relative to the center of the twin plane. This corresponds to the fact that 
the flux line is "free" and therefore completely delocalized in the direction parallel to the 
twin plane (x), while it is localized by the pinning potential in the direction transverse to 
the twin plane [y). The corresponding transverse localization length ly{T) can be deflned 
as 

/ + 00 
dyy'\My)\^ (2.7) 

-oo 

where the wavefunction is assumed to be normalized. As shown in pi|, one finds ly{T) ^ 
bo[l + 0{T/T*)], for T « T* . At high temperature thermal wandering is important and 
the localization length can become larger than 6o: with ly{T) ^ , _ ~ bo{T/T*)'^, 

for T>>T*. 

When many vortices are present, the repulsive intervortex interaction tends to confine 



each fiux line to a "cage" provided by the surrounding vortices in a triangular lattice [|2^ . 
For fields below the "filling" field Bf ^ (pg/cP, defined as the field where the fiux lines 
fill the twin planes, forming a triangular lattice of spacing \/3d/2 ||S^, the additional 
confining potential provided by the repulsive interaction does not change qualitatively the 
fiuxon states in the direction transverse to the twins, since in this direction fiux lines are 
also localized by the pinning potential. Interactions with neighbors do, however, change 



qualitatively the behavior along the x direction. Following Ref. [29|, a simple description 



of the role of interactions can be obtained by considering a representative fluxon localized 
near a single twin plane centered at y = and subject to the additional confining potential 
provided by the surrounding vortices. If the position of all the other vortices is assumed 
to be fixed, the confining potential can be approximated by a one-body effective potential 

Veffiriiz)). 

K//(ri(^)) = ^5^n|ri(^)-r;|), (2.8) 

where V{r) is the pair potential given in Eq. ( ^.5| ) and ri(z) denotes the position of 
the representative fiuxon. The sum is over all the other vortices that are fixed at their 
equilibrium positions r^, corresponding to the sites of the triangular lattice. The free 
energy of the representative fluxon is then given by, 



^1 



e// 



dz 



£1 
2 



dri{z) 



dz 



+ V,ff{r,{z))+V,{Mz)\) 



(2.9) 



where Vi(||/|) is the single- twin pinning potential discussed earlier. If we expand the 
effective potential Ve//(ri) about its minimum at ri = 0, we find 



K//(ri)^K//(0) + -Cr2 



(2.10) 



where, neglecting logarithmic corrections and constants of order unity, 

2eo 



C 



a, 



2 ' 



for Xab » Go, where the pair interaction V{r) is logarithmic {Ko{x) 

and 

C~ 2eo / 7rAab ^_^^/^^^ 



Kb^ 2ao 



(2.11) 
Inx, for X << 1), 

(2.12) 



for Xab « do, where the pair interaction decreases exponentially with distance. Again, 



following Ref. [^ , in the limit L — > oo the partition function of this representative fiuxon 
is written in terms of the ground state eigenfunction and eigenvalue of the "Hamiltonian" 
operator of a fictitious quantum mechanical particle. Dropping the constant term in Eq. 
( p.lOp , the corresponding "Schrodinger" equation is given by 



-^V? + ^Cr? + yi(|yi|)' 



^o(ri: 



^o^o(ri; 



(2.13) 



The X and y degrees of freedom are decoupled and the "Schrodinger" equation (|2.13|) can 
be separated into two one dimensional equations (to simplify the notation, we drop the 
subscript 1 on the location of the representative fluxon), 



2e~idx2 + 2 



g^{x) = E^g^{x), (2.14) 



and 

- §^^ + \cy^ + V^{\y\)\gy{y) = Eygy{y), (2.15) 

with ^o(r) = gx{x)gy{y) and Eq = Ej^ + Ey. In the x direction the vortex line is described 
by the ground state of a one- dimensional harmonic oscillator of frequency ujq = \/C jlx ~ 
(l/ao)A/2eo/ei. The ground state energy is E^ = Tloq and the corresponding eigenfunction 
is 

j{x) = —^ ^g-(x/2x*)2 ^2.16) 

^ ^ (V^a;*)i/2 ' ^ ^ 



where x* = (T2/2eiC)^/^ is the characteristic length scale for vortex fluctuations along the 
X direction. In the absence of the twin plane, the ground state in the y direction is also that 
of a harmonic oscillator of frequency wq and the vortex is confined by interactions within 
a region of radius r* = a/x*^ -f y*2, with y* = a;*, centered at its equilibrium position, 
r = 0. The presence of the twin boundary modifies the potential in the y direction, leading 
to an additional square well near the center of the harmonic potential, as in Eq. ( |2.15|) . 
The range of the wavefunction g{%j) controls the localization length Ix. in the direction 
transverse to the twin plane. This is determined by the interplay of the length scale 
y*{T) for harmonic fiuctuations and the localization length ly{T) defined in Eq. ( p.7| ) 
associated with the pinning potential. These two length scales are sketched in Fig. 5 
as functions of temperature. For T « T* , ly ^ Bq. If the temperature is so low that 
y* < ly ^ bo, the range of the wavefunction g{y) is controlled by interactions and l± ~ y*. 
The characteristic temperature Tj;i where y* = bo is given by Tj;i = (bo/ao) \f2eofUoT* « 
T*, as shown in Fig. 5. For T >> T*, /j/(T) ^ bo(T/T*)'^ grows more quickly than 
y* with temperature. There is therefore a second crossover temperature Tx2-, as shown 
schematically in Fig. 5. For T > Tj;2, ly > y* and the vortex line is confined only by 
the harmonic well from intervortex interactions. A lower bound for Tx2 can be obtained 
from y*(T,2) ^ ly ^ bo{T^2/T*f. with the result T,2 = [{ao/bo)^U^j2roY/^T* > T\ 
For the parameters of interest here Tj:2 is somewhat smaller than the clean lattice melting 
temperature Tm, defined by y*(Tm) ~ clOo, with cl ~ 0.15 — 0.3 the Lindeman constant, 
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and is comparable to the isolated vortex depinning temperature, defined by ly{Tdp) ~ d. 
In this paper we only consider the situation where all vortices are pinned in the ground 
state and T « min(Tdp, Tm)- We therefore restrict ourselves to T < Tj:2 and we then find 
l^ ~ y* for T < T^i and l±_ ~ ly for Tj^i < T < Tr^2- In short, the ground state of a single 
vortex confined by the pinning potential of a twin plane along the y direction and by the 
isotropic "cage" provided by the repulsive interaction with the other vortices is localized 
in all directions, with localization lengths l± ~ ly in the direction transverse to the twin 
(for Txi < T < Tx2) and /y ~ x* in the direction parallel to the twin plane. The total 
binding free energy renormalized by interactions is of order Ur{T) k, U{T) — Tujq. For 
T « T^2 the harmonic oscillator zero point energy is always negligible compared to the 
pinning energy U{T) and Ur{T) ^ U{T). 

In the presence of a family of parallel twin planes the flux line can "tunnel" between 
different localized states [p9| , p!3[] . In this paper we are interested in studying the response 
of a flux array pinned by a family of parallel twin planes to a Lorentz force normal to 
the twin planes for B « Bj. In the ground state the flux lines are all localized on 
the attractive twin planes. The transverse driving force promotes motion of the vortices 
between different twin planes, corresponding to "tunneling" between different localized 
states along the y direction, while the repulsive interaction confines the vortices in the 
direction parallel to the twin planes. At low temperature the flux lines will move along 
the direction of the driving force {y direction) within one-dimensional channels of width 
~ 2x* . Using elementary quantum mechanics it can be shown [p] , p9| that the rate of 
tunneling between localized states on different twin planes separated by a distance dij is 
tij ~ 2U{T)e-^^i/'^, with Eij = ^2iiU{T)dij. The energy E,j is the energy of a "kink" 
configuration shown in Fig. 3b, connecting two pins at a distance dij. 

To study the low-lying excitations from this ground state arising from thermal fluctu- 
ations one needs to sum over vortex trajectories by evaluating appropriate path integrals. 



As discussed in [jT^ and [^ , these conflguration sums closely resemble the imaginary time 
path integral formulation of quantum mechanics of two-dimensional particles in a static 
random potential Voiy)- Many relevant results regarding the statistical mechanics of flux 
lines can then be obtained from elementary quantum mechanics. 

The dynamics of flux lines driven by a Lorentz force transverse to the twin planes can 



then be described by a tight-binding model for one- dimensional hosons [|T^. The lattice 
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sites in the model are defined by the M positions {Yi} of the twin planes and the tight 
binding Hamiltonian governing the dynamics in each one-dimensional channel is given by 

if = -[// + U{T)]Y,a\a, + E^^^(4a. + 4«^) + y E4«^(4«^ - 1)- (2-17) 

i i^j i 

Here n ~ (po{H — ifci)/47r is the chemical potential which fixes the flixx line density, al 
and Qi are boson creation and annihilation operators at site Yi, tij is a tunneling matrix 
element connecting localized states i and j and Vq represents a typical energy cost for 
double occupancy of a site of the one-dimensional tight-binding lattice. As flux lines move 
in the transverse direction along the one-dimensional channels, the repulsive intervortex 
interaction provides an energy cost for an additional flux line occupying an already filled 
twin. The corresponding on-site repulsion Vq can be estimated as 

^0 ~ V{bo) - V{d) 

(2.18) 
^ 2eo[H>^ab/U) - Ko{d/Xab)], 



where V{r) is the pair interaction given in Eq. (|2.5| ) and we assumed bo ~ ^ab- li d >> Aa6, 
we find Vo ~ 2eo ln(Aa&/^a6), while for d « Aa6, we obtain Vq ~ 2eo ln((i/^a&)- This 
estimate assumes that the various one-dimensional channels are completely decoupled. 

The first two terms of the tight-binding Hamiltonian determine a noninteracting den- 
sity of states g{e) (here e is an energy per unit length and g{e) has units of 1/energy), such 
that J\f{e) = f_^ g{e')de' is the number of localized states per unit length with energy less 
than e. Note that g[e) is normalized so that A/'(-|-oo) = 1/d. Even if the pinning sites 
are all identical in size and well depth, dispersion of energy levels arises because vortices 
can tunnel between nearby twin planes. The width 7 of the impurity band should then be 
of order 7 k, t{d), where t{d) is the tunneling matrix element evaluated at a typical twin 
spacing d. Interactions will further broaden the band and one can estimate, 

7^max{t(rf),"l/o}- (2.19) 

This bandwidth is practically always dominated by Vq. In particular for the case d << Xab 
one finds 7 ^ Vq ~ 2eo ln{d/^ab) for all temperatures T < T*d/bo. If the localized states 
are filled up to a chemical potential fi such that about half of the twins are occupied by at 
least one vortex, we can approximate the density of states g{fx) corresponding to the most 
weakly bound fiux lines with energy e ~ /i as g{n) ~ l/(i7, i.e., 

dgi^i)^mm{l/t{d),l/Vo}. (2.20) 

From our discussion of the bandwidth 7 we find that the second term generally dominates 
for all temperatures of interest. Then g{ij,) is approximately temperature independent, 
with dg{fx) ~ I/Vq. 

12 



3. Vortex Dynamics at Low Temperatures 

We consider vortex transport in tlie presence of a driving current J ± H parallel to 
the twin planes, i.e., J = —Jit. The applied current exerts a Lorentz force per unit length 
on the vortices (see Fig. la), 

fL = -zxJ = yh, (3.1) 

c 

with Jl = (poJ/c, and drives the vortices to move in the direction transverse to the twin 

planes, leading to an additional term, 

SFi = ~fL [ yi{z)dz, (3.2) 

Jo 

in the single- vortex free energy, Eq. (|2.1|) . In the context of the analogy with boson 
quantum mechanics, this term represents a fictitious "electric field" E = -z x J = yJ/c 
acting on particles with "charge" cpQ. The correspondence between the problem of carrier 
dynamics in disordered semiconductors and vortex dynamics in the presence of correlated 
linear or planar disorder is summarized in Table 2. 

Up to numerical constants and logarithmic corrections, the critical current at low 
temperatures can be obtained by equating the Lorentz force to Uo/bo, with the result 
Jc{0) ~ cUo/(pobo [|T3| . Thermal fluctuations renormalize the critical current and one can 
estimate J^T) ^ cU {T) / 4>qI ^{T) . For T » T*, J^{T) ^ J^{'d){T* /Ty. The crossover 
temperature T* is itself a function of temperature. We can define the temperature Ti above 
which the entropy from fiux-line wandering is important in renormalizing the binding free 
energy U{T) by the self- consistency relation T*(Ti) = Ti[|21|]. At low temperature the 
interaction between a vortex line and a twin plane is always attractive and Uq ~ a^toT, 
where ctf, < 1 is a dimensionless parameter related to the barrier transparency and r = 
1 — T/Tcpl]. Using a mean field parametrization of the critical fields, we find T*{T)/Tc = 



i^r, where k = Kb/U and Gi = {\l/2\l^){TjHl^eaw? is the Ginzburg number, 
with Hco the thermodynamic critical field at T = and ^ah = ?ab^~^^^- Using k k, 10^, 
Gi ^ 10-2 and af, ^ 0.1 g for YBCO, we find Ti ^ 0.77Te. 

At low temperatures and fields well below Bf, vortex dynamics is determined by the 
competition between pinning by the one- dimensional array of twin planes and thermal 
fiuctuations of the vortices. In analogy with the case of columnar pins that was discussed 
in detail in [jl^ , the boson mapping reduces single vortex dynamics to a problem of hop- 
ping conductivity of localized particles in one dimension. The current density in the usual 
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hopping conductivity problem corresponds to the vortex velocity (i.e., voltage) and the 
electrical conductivity maps onto the resistivity from vortex motion (see Table 2). The 
low temperature dynamics of vortices driven transverse to an array of parallel twin planes 
presents the same rich variety of hopping conductivity phenomena that occur in semicon- 
ductors, as pointed out by Nelson and Vinokur for the case of columnar pins fl^ . What is 
new here is that vortex dynamics maps onto the problem of hopping conductivity in one 
dimension. In this reduced dimensionality rare events, such as large regions voids of twin 
planes, can dominate the transport at low currents leading to new mesoscopic phenomena, 
as discussed in Section 5. 

Here we are interested in the low temperature regime where transport is dominated 
by single-vortex dynamics. In this case the dominant contribution to dissipation can be 
described in terms of the low-lying excitations from the ground state that correspond to 
thermally activated jumps of vortex lines over the relevant pinning energy barriers. The 
resistivity takes the form given in Eq. ( |1 . 1| ) . In the following we determine the barrier 
heights U{L, J) corresponding to various transport regimes and the boundaries between the 
various regimes in the (L, J) plane. In samples of finite thickness L in the field direction the 
typical pinning energy barriers W(L) grow with L but are independent of current, yielding 
a linear resistivity. In thick samples there is a nonlinear resistivity associated with barriers 
U{J) that grow at low currents. We assume that in the ground state all the flux lines are 
localized on twin planes. We then study the low lying excitations from the ground state 
that can be nucleated by a flnite temperature T or by a driving current J. The largest 
contribution to the resistivity from each class of excitations is assumed to be inversely 
proportional to the shortest time for the nucleation of a given excitation. The latter is 
determined by the typical energy barrier U for the formation of the low-lying excitations, 
which is identifled with the saddle point in the single-flux line free energy. The discussion 



in this section follows closely that of Ref. [|T3[ , where the corresponding results for votices 
pinned by columnar defects were obtained. 

Linear response 

Consider a fluctuation that extends a length z along the twin and a distance y in 
the direction of the Lorentz force. The free energy of this fluctuation relative to the case 
/l = is 

SF{y,z)^ei^ + Uz-fLyz. (3.3) 
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Optimization of Eq. ( |3.3|) with repsect to z for /l = yields the shape of the optimal low 
temperatute fluctuation, 

y ~ ^JWJTx z. (3.4) 

At low currents in samples of very small thickness L there is a linear resistivity due to the 
flow of rigid flux-line segments of length L and typical transverse width y^/ ~ \/U{T)/eiL 
obtained by letting z ~ L in Eq. (|3.3|). The corresponding saddle point energy is lArf{L) ~ 
U(T)L, resulting in a linear "rigid-flow" resistivity, 

PrfiL)^poe-^(^^^/^. (3.5) 

At larger currents the contribution from the Lorentz force to the single-line free energy 
( ^73|) becomes comparable to the typical energy barrier Urf- When fLVr/L > Urj or 
J > Jl = c\^eiU/{(poL), the response becomes nonlinear. In the thermodynamic limit 
Jl —>■ and the IV characteristic is nonlinear at all currents. The characteristic current 
Jl is also conveniently expressed in terms of the energy of a kink conflguration connecting 
neighboring pins separated by the distance d (see Fig. 3a). The typical thickness Wk 
of a kink along the z direction is obtained from Eq. ( |3.4]) for y r^ d, with the result 
Wk = d\Jl\fU . The kink energy is the corresponding saddle-point free energy, E^. = 
WkU = \JiiU{T)d. The energy barrier associated with the "rigid-flow" resistivity can 
then be written as Urf = Ek{L/wk) and the current scale for nonlinear transport is Jl = 
cEk/{4>oLd). 

The line J = Jl{L) deflnes the boundary in the (L, J) plane that separates the regions 
of linear (J < Jl{L)) and nonlinear (J > Jl{L)) response (see Fig. 2). The details of 
the (L, J) phase diagram are controlled by the dimensionless parameter a = g{p)dU{T). 
Typical phase diagrams for a < E^/T are shown in Fig. 2. The rigid flow mechanism 
dominates the linear resistivity only in very thin samples. When Wk < L < Li, where 
Li = Ek/'j is the length below which dispersion from tunneling and interactions can be 
neglected, transport occurs via the hopping of vortices between nearest neighbor (nn) 
pinning sites. This region of the phase diagram is only present if Li > d\/t\l\J , or 7 < t/. 
For L > L\ dispersion is always important and the relevant excitations are superkinks 
(Fig. 3c), which correspond to the tunneling of vortices between remote pinning sites 
analogue to Mott's electronic conductivity in disordered semiconductors. For J > Jl{L) 
the resistivity is nonlinear. At large currents the typical transverse displacement is smaller 
than the average spacing d between twin planes and transport is dominated by "half-loop" 
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excitations (Fig. 3a), characterized by an energy barrier that grows hnearly as the current 
decreases, Uhi ~ 1/ J . FinaUy, at the smaUer current flux motion takes place via VRH, 
characterized by a diverging energy barrier, IAvrh ~ 1/J^'^. We now discuss in more 
detail the origin of the various contributions summarized in Table 1 and the estimate of 
the energy barriers. 

We flrst consider the linear portion of the phase diagram (J < J^) in samples of 
increasing thickness L. When the typical transverse width yrf ~ \JU /ei)L of a rigidly 
flowing flux segment becomes comparable to d, transport occurs via nucleation of double 
kink conflgurations (Fig. 3b) of energy ~ ^E^. The double kink then separates to z = ±oo. 



resulting in the hopping of vortices from one pin to a neighboring one. As discussed in |13], 
this transport mechanism will dominate only if the sample is so thin that the width 7 of 
the impurity band arising from tunneling and interactions is negligible (L < Li = E^/'j). 
In this case flux motion will occur via hopping between nearest neighbor pins, resulting 
in a linear resistivity Pnnh ~ exp(— ai?fc/T), with a a numerical constant. In extremely 
thin samples this transport mechanism will ultimately be suppressed. In fact for L < Wk 
transport via the flow of rigid flux segments described above is energetically favorable over 
nn hopping. As a result, a necessary condition for observing a linear nearest neighbor 
hopping resistivity Pnnh is Li > Wk, or ■j < U. If we estimate the density of states as 
dg{fx) ~ 1/7, the condition ■j < U requires a = g{fx)dU > 1. The Lorentz force term in 
Eq. ( p.3|) will modify the kink energy and thickness. By optimizing Eq. ( p.3|) with respect 
to z for y ^ d and /^ 7^ 0, we flnd that a flnite current increases the thickness of a typical 
kink, according to Wk{J) = Wk{^ — J/Ji)~^/^, where Ji = cU/{(j)od). This result only 
applies for J < Ji. At higher currents simple nn hopping cannot occur. 

In thick samples (L > Li) the dispersion of energies between different pinning sites 
makes motion by nearest neighbor hopping energetically unfavorable (the energy barrier 
diverges with the sample thickness L). Tunneling occurs instead via the formation of 
"superkinks" (Fig. 3c) that throw a vortex segment onto a spatially remote pin connecting 
states which optimize the tunneling probability. The free energy of a superkink excitation 



shown in Fig. 3c relative to the case /^ = is then [|13| , |27 



5F,k ^ 2Ek (y/d) + Aez - -hyz. (3.6) 

We assume all states up to a chemical potential p are fllled. The states available to a weakly 
bound flux line about to hop a distance y are those within an energy Ae determined by 
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requiring that there is at least one locahzed state within a region (y, Ae) of configura- 
tion space, i.e., g{ii)yAe ~ 1. The shape of the most important superkink excitations is 
obtained by minimizing Eq. (|3.6|) for /l = and is given by 



In finite thickness samples the saddle point free energy corresponding to superkink fluc- 
tuations of width given by Eq. ( |3.7| ) for z ^ L yields a linear Mott resistivity, given 

by 

PMottiL) ^ poe-^"^"^/^-"^"'", (3.8) 

with a = U{T)g{fx)d. 

Nonlinear response 

In the nonlinear regime (J > Jl{L)) the contribution to the free energy from the 
Lorentz force cannot be neglected when estimating the energy of the dominant excitations. 
For Ji < J < Jc, with Ji = cU{T)/(pod, flux motion occurs via thermally activated "half- 



loop" conflgurations identical to those discussed in |jT^ for the case of columnar pins. The 
length and width of an unbound line segment for the lowest-lying half-loop excitations are 
obtained by minimizing the free energy (|3.3|) for fL 7^ 0, with the result Zhi ~ (t/^oei)^/^//^ 
and Phi ~ {Uo/ei)^^'^Zhi ~ U{T)/ Jl, respectively. The saddle point energy of a half-loop 



excitation is Uhi ~ y'eiU^{T)/ Jl, yielding a nonlinear resistivity. 

Phi ~ poexp[~{Ek/T){Ji/J)]. (3.9) 

In the context of the mapping of flux-line dynamics onto the problem of hopping conduc- 
tivity, the nucleation of half loops corresponds to tunneling of a carrier from a localized 
state directly into conduction band, as shown in Fig. 6. 

For J < Ji the size of the transverse displacement of the liberated vortex segment 
exceeds the average distance d between twin planes and transport occurs via variable range 
hopping (VRH) which generalizes the Mott mechanism to the nonlinear case. Again a flux 
line hops to a state within a region (y, Ae) of phase space, with g{p)yAe ~ 1. The size of 
the most important excitations is determined by minimizing Eq. ( |3.6| ) with f^ 7^ 0, with 
the result pvRH ~ (fi'(/^)/L)~^^^ cind zvrh ~ E^/dfL- One then obtains a non-Ohmic 
VRH behavior with 

PVRH ~ poe-"-«^(^)/'^ = poexp[-(i?fc/T)(Jo/J)i/2], (3.10) 
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with Jq = Ji/a. The crossover to the hnear Mott resistivity resistivity takes place when 
zvRH ~ i^, or J ~ Jl, consistent with the result obtained above when discussing half- loop 
excitations. The VRH contribution to the resistivity dominates that from half loop only 
if if UvRH <Uhh or J < J2 = Jia. 

The above results are summarized in Table 1. The corresponding phase diagrams 
are shown in Figs. 2 for a < E^/T. There are three relevant current scales, Jq = Ji/a, 
Ji = cU/(pod and J2 = Jict, all much smaller than the pair breaking current Jph = 
4ceo/(3-\/3(/)o^a6)- For a > 1 one can have Li > Wk and there is a region of the phase 
diagram where transport occurs via nn hopping (Fig. 2a). For a < 1 nn hopping can 
occur only if the chemical potential /x falls in the tails of the impurity band, so that 
g{y)d < 1/7. If n falls well within the impurity band, so that g{fJ.)d ~ I/7, then a < 1 
requires •y > U and nn hopping is always suppressed in this case. The Mott and the rigid 
flow regimes are separated by a horizontal line above which Umou < Urf- Similarly, the 
condition Uvrh = Uhi yields the vertical line separating the VRH and half loop regions. 

Collective effects 

At very low currents and in thick samples collective effects are always important 
and flux motion takes place via the creep of vortex bundles, rather than single vortices. 
The region where collective effects dominate is shown schematically in Figs. 2 and 4. It 
corresponds to the upper left portion of the (L, J) phase diagram. As discussed in Ref. 
0] , the crossover from single vortex creep to creep of vortex bundles occurs when 

Lz = ao, (3.11) 

where L^ is the size of a typical single-vortex fluctuations along the z direction and ao 
the intervortex spacing. The condition ( |3.11|) is simply obtained by equating the tilt 
energy Euu of a single disorted vortex to the elastic energy Eint of interaction of with 
its neighbors. The elastic energy associated with displacing a length L^ of vortex line at 
an average distance ao from its neighbors a distance u out of its equilibrium position in 
the xy plane is Ei^t ~ cggw^L^, where cee ~ ^qcLq, and grows with L^. In contrast, the 
corresponding single- vortex tilt energy, Emt ~ eo{u/Lz)'^Lz, decreases as L^ increases. 
Consequently when the longitudinal size of the typical fluctuation is sufficiently large, or 
L^ > ao, then Ei^t > Etm and collective effects are important. 

In the VRH regime the relevant length scale is the width Wsk of a superkink excitation, 
shown schematically in Fig. 3c, where Wsk ~ yhjUysk, with ysk the typical size of a 
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superkink in the direction of flux motion. Collective effects dominate when Wgk > flo- Iii 



the linear regime (J < Jl) Hsk is given by Eq. ( p.7|) with z ~ L, or Wsk ~ V^kL, with 
Wk = Ek/U the width of a kink (see Fig. 3d). Dissipation is then dominated by creep 
of vortex bundles for L > Lb = a^/wk- In the nonlinear regime (J < Jl) the size of the 
superkinks grows with decreasing current and Wsk ~ Wk{Jo/ JY^"^- Transport is always 
dominated by collective effects at sufficiently low currents, i.e., for J < Jb = Joiwk/ao)'^. 
The crossover from single- vortex creep to creep of vortex bundles is marked by the dashed 
lines L = Lb and J = J^ in Figs. 2 and 4. For a > 1 this crossover takes place well into 
the VRH region, as shown in Fig. 2a, provided Lb > Li and Jb < Ji, which corresponds to 
B < {Bf/a){U/2ei) (here and below we assume the chemical potential falls in the middle 
of the impurity band and dg{p.) ~ 1/7)- For a < 1 (Fig. 2b) this crossover occurs within 
the VRH region provided Lb > Wk/ct and Jb < J2, or S < Bfa{U/2ei) 

4. Transport in the Presence of Tilt 

We now consider another transport geometry investigated in some of the experiments 
by Kwok et al [Q. Here the external field H is tilted at an angle 9 away from the c axis 
and out of the twin planes (see Fig. lb). The transport current is still applied along 
the twin planes, which contain the c axis, J = — xJ, and the resulting Lorentz force, 
fi = {(po/c)[zcos9 + ysin^] x J, has components both normal to the twin planes and 
along the c axis. Only the y component of the Lorentz force is effective at driving flux 
motion normal to the twin planes and therefore determines the voltage in the direction 
of the applied current. The experiments by Kwok et al. Q have been mostly carried out 
at high fields for flux arrays in a liquid state, in a regime where intervortex interactions 
are believed to be important. Here in contrast we neglect intervortex interactions and 
investigate the dependence of transverse transport on tilt angle in the regime where single- 
line dynamics dominate. Even though our result are therefore not directly relevant to the 
experiments by the Argonne group, the strong angular dependence that we predict for the 
resistivity is qualitatively similar to that reported in the experiments. 

The free energy of a fluctuation that extends a length z along the twin and a distance 
y in the direction of average motion is obtained by adding the tilt energy Eq. (|3.3| ), with 
the result, 

dF{y,z,9) ^i,y- + Uz- U cos 9y z - ^H^y, (4.1) 

z 4n 
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where H±^ = if sin 6* is the component of the field along the y direction and /^ = (poJ/c, 



as in the preceeding sections. Optimizing Eq. (|4.1|) with respect to z for /^ = 0, we find 
that the shape of the optimal low temperature fiuctuation is still given by Eq. ( p^) and 



does not depend on the angle 6. As shown by Hwa et al. ||2^, the energy Ek{0) of a kink 
fiuctuation in the presence of tilt, corresponding to the saddle point of the free energy 
( ^?T|) with /i = for y ~ (i and z ~ d^/U/ii, is reduced compared to its value for ^ == 0, 
according to, 

^''. , (4.2) 

eJi ""^^ 



sin 9c ■ 

Here we have introduced a critical angle 9c defined by sin 9c = H/ He, with He = 4TTEk/(/)od. 
For 9 > 9c, the kink energy becomes negative and kinks proliferate, as discussed in |27[ j. 
We now consider the angular dependence of transport for 9 < 9c- 

As in the case ^ = 0, at high enough currents fiux motion will occur via the nucleation 
of half loops. By identifying the typical energy barrier Uhi{9) for a half loop excitation in 
the presence of tilt with the saddle point energy found by minimizing Eq. ( [4.1| ) for f^ 7^ 0, 
we obtain, 

".,W = ^(l-«^), (4.3) 

COS 9 V sm 9c / 

with a a numerical constant of order one and Uhi = Ek{J/Ji) the half loop energy barrier 

for ^ = 0. The energy for nucleating a half loop excitation is reduced by the tilt. The 

angular dependence of the resulting fiux-fiow resistivity is very strong, since the angle 

appears in the argument of the exponential. Flux motion will occur via half loop excitation 

provided the typical transverse size of the half loop does not exceed the average distance 

between twin planes. This imposes a lower bound on the values of the current where 

half loop excitation dominates transport, given by J > Ji(l + a' sin ^/ sin ^c) 5 with a' a 

numerical constant of order one. Tilt decreases the range of currents where half-loops 

dominate. 

At lower currents transport will take place via VRH. The angular dependence simply 

replaces the kink energy Ek by the smaller kink energy Ek [9) in the presence of tilt given 

in Eq. (|4.2|) . Carrying then through the standard VRH argument described in the previous 

section, one obtains a nonlinear angle-dependent resistivity given by Eq. ( |1 . Ij ) , with 

WvK«(9) = £.(Ji/«/)'/^i^^^^£^. (4.4) 

vcosy 
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Again the corresponding resistivity is a rapidly varying function of angle, as observed in 
experiments. On the other hand, a simple estimate using typical parameters for YBCO 
gives a very small value for the critical angle, sin^^ ~ O.lHd/H. The transport experi- 
ments probe, however, linear transport in the flux-liquid phase, where collective effects in 
the flux-line dynamics are important. The present dimensional analysis is useful in that it 
shows that even in the regime of single-line dynamics, the presence of twin planes naturally 
introduces a very sharp dependence of the resistivity on tilt angle. 

5. Rare Fluctuations 



The results described in the Section 3 are qualitatively similar to those discussed in ^3 
for the case of flux arrays in the presence of columnar pins. The most important difference 
for samples with parallel arrays of twin planes is that due to the one-dimensional nature 
of vortex transport at low temperature, a new regime can arise at low current, where flux- 
line dynamics is dominated by rare fluctuations in the spatial distribution of twin planes. 
The vortex line can encounter a rare region where no favorable twins are available at the 
distance of the optimal jump. The vortex will then remain trapped in this region for a 
long time and the resistivity can be greatly suppressed. Rare fluctuations can also occur 
in samples with columnar pins, but in that case because of the two-dimensional nature of 
the problem, they will dominate transport and suppress the resistivity only at extremely 
small fields, when the number of rare regions exceeds the number of vortices. 

At a given temperature and for applied currents below J^, a vortex can jump from 
one twin plane to another at a distance y only if the energy difference per unit length 
between the initial and final configuration is within a range Ae ~ Eky/Ld. A trap is then 
a region of configuration space (y, e) void of localized states within a spatial distance y and 
an energy band Ae around the initial vortex state. A vortex that has entered such a trap 
or "break" will remain in the trap for a time t^^ ~ to 6xp(22///^), where l± is the transverse 
localization length and to is a microscopic time scale. The probability of finding such a 
break is given by a Poisson distribution, P{y) ~ Po{y)exp[—Ag{fx)yAe], where Po{y) is 
the concentration of localized states in the energy band Ae, Po{y) ~ 2Ag{^)Ae and A ~ 1 
is a numerical constant. The mean waiting time between jumps is given by 



/•OO 

/ dyP{y)toe^y/'-^^l (5.1) 
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For L >> L* = awk{T/Ek)'^, the integral can be evaluated at the saddle point, corre- 
sponding to the situation where the mean waiting time is controlled by "optimal breaks" 
of transverse width yf ^ l±L/L* , with the result, 

t;:~to/W^e^/^*. (5.2) 

The optimal breaks are those that correspond to the longest trapping time and will there- 
fore be most effective at preventing flux motion and dissipation. The inverse of the trapping 
or waiting time determines the characteristic rate of jumps, i.e., the velocity. The vortex 
velocity corresponding to the optimal hopping rate of Eq. ( |5.2| ) yields a linear resistivity 
in finite-thickness samples, given by 

Pw~P077^e-("-^°"/'^)', (5.3) 

l^Mott 

where Umom is given in Table 1 and we have used L/ L* = {Umom/TY . 

For currents above J^, the typical energy per unit length available to a flux line 
for jumping a distance y is Ae ^ JlU- The corresponding nonlinear contribution to the 
resistivity from traps of extent (y, Ae) in configuration space is again proportional to the 
inverse of the average waiting time defined in Eq. ( |5.1j ). Again, for l\g{^)fL << 1 or 
J « J* = aJi{Ek/T)'^, the integral can be evaluated at the saddle point, corresponding 
to an optimal break width y^ ^ [g{fJ')i±fL]~^j with the result, 

Pb ~ pori^exp[~{UvRH/Tf]. (5.4) 



It is clear by comparing Eqs. (|3.9| ) and (|3.10|) to Eqs. (|5.3|) and ( p.4|) , respectively, that the 



contribution to the resistivity from tunneling a la Mott (both in the linear and nonlinear 
regimes) would always dominate that from hopping between rare optimal traps if both 
mechanisms of transport can occur. On the other hand, in one dimension if the sample is 
wide enough in the direction of flux-line motion to contain optimal traps, tunneling a la 
Mott simply cannot take place because flux lines cannot get around the traps. These rare 
traps with large waiting times will then control the transport. If W is the sample width 
in the y direction, the condition for having optimal traps of width y^*^ is P(y;*^)VF > 1. 
Optimal traps will therefore be present only if J > J^ = J* /\n{2W/l±) for J > J^ and if 
L<L^ = L* \n{2W/l±) for J < Jl- 

These are, however, only necessary conditions for the sample to contain many optimal 
breaks. They do not guarantee that these breaks will dominate transport. A flux line can in 

22 



fact escape a break by nucleating a half- loop excitation or, in the language of semiconductor 
transport, by tunneling directly from a localized state into conduction band (see Fig. 6). 
This will occur if the transverse size of a typical half loop exceeds the size of the optimal 
break, i.e., if yhi > yl, or a < E^/T. In the context of the analogy with boson quantum 
mechanics this condition translates into the requirement that the spatial distance between 
the occupied localized state and the conduction band edge in the presence of the applied 
current is shorter than the size of the trap (see Fig. 6). If the flux line can escape the trap 
by half-loop nucleation, breaks will never dominate transport and their only effect will 
be that of possibly suppressing VRH in a region of the phase diagram and extending to 
lower currents the region where transport occurs via half-loop nucleation. This can occur 
if J-w is smaller than the scale setting the high current boundary of the VRH region, or 
Ju, < min(Ji, J2). For instance if Jyj < J2, or \n{2W/l^) > (Ek/aT)^, with a < 1, rare 
fluctuations will modify the phase diagram of Fig. 2b by pushing the high current boundary 
of the VRH region down to Jw Similar considerations apply to the linear response. On 
the other hand, if yhi < yl, or a > E^/T, there will be a portion of the {L,J) phase 
diagram where breaks dominate transport, as shown in Fig. 4. 

For YBCO, we estimate E^ ~ lKA~^d. Assuming a ~ U/'j, the condition a > E^/T 
can only be satisfied at low fields {B < IKG for d ~ 200A). The sample will contain 
optimal breaks if W > 30Aexp(J*/J), with J* ~ 4 x lO^Amp/cm^ at 80K. 

If the sample is too short to contain optimal breaks, i.e., WP{y ~ yf ^) < 1, the 
dynamics is controlled by the trap with the longest waiting time, t{yf) ~ exp(yj//^), with 
yf determined by the condition WP{yf) ~ 1. The corresponding resistivity is proportional 
to this smallest hopping rate, 

Pw^Poe-"^'/'^. (5.5) 

In this case the relevant physical quantity is the logarithm of the resistivity. 



Wv 



21^ T / „ aV2 



1/2 



HPw/Po) = -yf/i^^-^{\^[-j^jj^[H'^w/iS) ]} . (5.6) 

The leading dependence of Eq. ( |5.6|) on current and temperature is the same as that 
of the VRH contribution. Equation (|5.6|) also contains, however, logarithmic terms that 
in sufficiently short samples will give a random spread of values of the resistivity from 



sample to sample. These effects have been discussed for semiconductors [jT^. In this case 



a more relevant physical quantity rather than the resistivity itself is the distribution of the 
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logarithms of the resistivity over different samples. The expression ( ^.61) determines the 
position of the maximum of this distribution. 
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Linear 


Nonlinear 


Urf = UL = Ek{L/wk) 


UHi = Ek{Ji/J) 


Unnh=Ekil-J/Jl)-'/^ 




V(Mott=EkiL/awky/^ 


UvRH = EkiJjajy/^ 



Table 1. Energy barriers determining the various contributions to the resistivity of Eq. 
(PI), Mith a ^U{T)g{^)d. 
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CARRIERS 


VORTICES 


m 


ei 


h 


ksT 


pn 


L 


single impurity level Ed{Ea) 


U{T) 


fi 


{cPo/4n){H-H,i) 


E 


c 


carrier velocity ~ current density 


vortex velocity ~ voltage 


conductivity a 


resistivity p 


conduction-band transport 


flux flow 


tunneling from impurity levels to conduction band 


half- loop 


VRH 


superkink 



Table 2. Correspondence between carrier dynamics in disordered semiconductors and 
vortex dynamics in the presence of correlated disorder. 
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Figure Captions 

Fig. 1. Geometry of the transport experiment corresponding to strong pinning by twin 
boundaries. In (a) the external field is aligned with the c axis and lies in the 
plane of the twins. In (b) the external field is tilted at an angle 9 out of the plane 
of the twin. In this case only the y component /l cos 9 of the Lorentz force is 
effective at driving flux motion transverse to the twins. 

Fig. 2. The (L, J) phase diagram for a = g{fx)dU ~ t//7 < E^/T. The curved phase 
boundaries between the Mott and VRH regimes and between the rigid flow and 
half loop regimes are determined by Jl = cEk/(podL. For currents below the 
characteristic current scale Ji the typical transverse size of a fluctuation in the 
nonlinear portion of the diagram exceeds the average distance d between pins. 
The corresponding length scale Wk = Ek/U is the width of a kink connecting 
pins at the distance d. The crossover from half-loop to VRH is determined by 
min(Ji, J2), with J2 = aJi. Figure 2a is for a > 1, corresponding to Li > Wk, 
where Li is the sample thickness above which level dispersion is important. In 
this case flux motion can take place via nn hopping for Wk < L < Li and J < Ji. 
Figure 2b is for a < 1, when Li < Wk and nn hopping is suppressed. The typical 
energy barriers determining the resistivity in the various regimes are given in 
Table 1 . The dashed lines in the upper left corner of the plane delimit the region 
where collective effects are important (see text). 

Fig. 3. Schematic representation of the various low-lying excitations discussed in the 
text: (a) half-loop excitation, (b) double-kink conflguration, with Wk = (iyfi/Z/, 
and (c) double-superkink conflguration required for VRH. 

Fig. 4. The (L, J) phase diagram for a = g{n)dU > Ek/T. In this case there is a region 
where the TAFF resistivity is controlled by rare regions, both above and below 
Jl- The width of this region is controlled by the sample size W in the direction 
of flux motion. 

Fig. 5. The localization lengths in the direction transverse to the twin planes: ly(T) 
is determined by the pinning potential of the twin and y*{T) is determined by 
intervortex interactions. The various temperature scales are discussed in the text. 

Fig. 6. Schematic skecth in conflguration space (e, y) illustrating that a half-loop ex- 
citation corresponds to tunneling of the flctitious quantum mechanical particle 
directly into conduction band. The straight line of slope —fL/U{T) is the con- 
duction band edge in the presence of the flctitious electric fleld due to the Lorentz 
force. A carrier occupying a localized state at y = near the center of the im- 
purity band, i.e., at an energy ~ U{T) below conduction band edge, is brought 
directly into conduction band by a hop of transverse size yhi ~ U{T)/ f^. 
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